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Abstract. We explore a cosmological model composed by a dark matter fluid interacting with 
a dark energy fluid. The interaction term has the non-linear Xp^Pe form, where p m and p e 
are the energy densities of the dark matter and dark energy, respectively. The parameters a 
and f3 in principle are not constraint to take any particular values. We perform an analytical 
study of the evolution equations, finding the fixed points and their stability properties in 
order to characterize suitable physical regions in the space of the dark matter and dark 
energy densities. The constants (A, a, (3) as well as (u; m ,-u; c ) of the EoS of dark matter and 
dark energy respectively were estimated using the cosmological observations of the type la 
supernovae data set and the Hubble parameter H{z) at different redshift. We found that 
the best fit to data is for a model with a phantom dark energy interacting with a warm dark 
matter, where the energy transfer comes from the dark energy to the dark matter and with 
an interacting term of the simple form "/9 m p c " . This result is consistent with stable solutions 
of the dynamical system analysis. 
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1 Introduction 

The existence of a dark component with an exotic equation of state, i. e., with a ratio 
w = p/p negative and close to —1, which drives an accelerated expansion, is consistent with 
data coming from Supernovae la (SNe la) [1], large scale structure formation (LSS) [2], cosmic 
microwave background radiation (CMB) [3], baryon acoustic oscillations (BAO) [4] and weak 
lensing [5]. 

Cosmic observations show that densities of dark energy (DE) and dark matter (DM) 
are of the same order today, despite their different decreasing rates. To solve this coincidence 
problem [6] an evolving dark energy field with a non-gravitational interaction with matter [7] 
is assumed (decay of dark energy to dark matter). In this case both dark fluids interact via 
an additional coupling term in the fluid equations. In the current research different forms of 
coupling have been considered. Most of them study the coupling between cold DM and DE. 
In general, the interactions investigated are particular cases the form X m H p m + X e H p e , where 
H is the Hubble parameter, and p m and p e are the dark matter and dark energy densities 
respectively. Nevertheless, non linear interactions of the form A Pm £ e were considered in [8]. 



A more plausible interaction is inspired by the situation of two types of fluids interacting, 
where the interaction is proportional to the product of the powers of the number density of 
both components. In this case the interaction rate goes to zero as one or both densities become 
zero, and increase with each of the densities. We then consider in this paper an interaction 
of the form Xp^Pe, where the parameter A has dimensions of [density a+l3 ~ 1 x time] -1 . This 
type of interaction was investigated in the framework of an holographic dark energy [9] and in 
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coupled quintessence [10], where the evolution of the energy densities of the dark interacting 
components was investigated for different values of the parameters a and j3. A cyclic scenario 
for the present situation p e 

^ Pm was found as a possible solution to the coincidence problem. 
The particular case a = (3 = 1 was studied in [11]. In this case if the energy is transferred 
from dark energy to dark matter (A > 0), and for a dark energy of a phantom type (w e < — 1) 
the energy densities also presents periodic orbits. A general result, which is independent of 
the interaction type, and only assumes that the energy is transferred from dark energy to 
dark matter, pointed out that stationary solutions of the ratio r = p m /p e require a phantom 
dark energy [12]. 

The aim of this paper is to study the time evolution of the dark sector densities when 
the above mentioned interaction is considered. This evolution is driven by a highly non-linear 
coupled differential equations when the parameters a and /3 are let free. We solve numerically 
these equations and also study analytically the stability of the fixed-points. Since has been 
argued that second law of thermodynamics and Le Chatelier's principle implies A > [13] 
we do not consider here the case of energy transferred from dark matter to dark energy. 

In most of the investigation which take into account this type of non lineal interaction, 
cold dark matter with zero pressure is taken to be interacting with the dark energy fluid. 
Nevertheless, current research has opened the possibility that a warm dark matter component 
fits better new results found at the level of galaxies and cluster of galaxies [14]. For a wide 
discussion on this matter see [15]. In this work we investigate how a warm dark matter 
modifies the behavior of an interacting dark sector and the constraints for its EoS (assuming 
a barotropic form) derived from astrophysical observations. 

Our paper is organized as follows. In section II we present the effective equations of 
state for the interacting dark components, assuming that the energy is transferred from dark 
energy to dark matter. We also extend the result found in [12] in order to find the condition 
for stationary points of the parameter r when the dark matter is assumed not to be dust. 
In section III we study the two coupled differential equations corresponding to continuity 
equations of both interacting fluids. We study analytically the fixed-points including their 
stability. In section IV, we obtain numerical solutions using a numerical method of adaptive 
step-size algorithm called Bulirsch-Stoer method. We explore the behavior of the fixed points 
varying the parameters a, f3 and the EoS of both dark sector components. 

We study the stationary solutions and their stability for the equations of evolution of 
the ratios between the DM and DE and the third fluid and DE. We present the constraints on 
the equations of state for DE and the hypothetical third fluid. In section IV we discuss our 
results. Section V is an Appendix in which the conditions for the stability of the stationary 
solutions, corresponding to the three different couplings, are discussed. 

2 General considerations of an interacting dark sector 

In the following we assume a flat FRW universe filled basically with the fluids of the dark 
sector. We consider a warm dark matter of density p m and a dark energy component described 
by the density p e . For simplicity we also assume that both fluids obey a barotropic EoS, so 
we have p m = w m p m , for the warm dark matter and p e = w e p e for dark energy. In what 
follows we restrict our model to the late time of cosmic evolution, which implies that the 
others components of the universe, like radiation and baryons are negligible. In this case the 
sourced Friedmann equation is given by 

3H 2 = Pm + p e , (2.1) 
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where 8irG = 1 has been adopted. We will assume that the dark matter component is 
interacting with the dark energy component, so their continuity equations take the form 



p m + 3(1 + w m )Hp m = +Q (2.2) 
p e + 3(1 + w e )Hp e = -Q, (2.3) 

where H = a/a is the Hubble parameter, and a(t) is the scale-factor. Here, an overdot 
indicates a time derivative. Q represents the interaction term, despite that we do not use a 
specific functional dependence at this stage, we will only assume that Q do not change its 
sign during the cosmic evolution. 

2.1 Effective EoS for interacting fluids 

Let us discuss briefly the behavior of the dark components in terms of an effective EoS driven 
by the interacting term. Rewriting the continuity equation (2.2) in the usual form 

p m + 3(1 + Wm E}} )H Pm = 0, (2.4) 

where w meff represents the effective EoS for the interacting dark matter, which is given by 

Wm tt = w m — • (2.5) 

m eff m 3Hpm V ) 

Note that the behavior of w meff can be quite different if w m is non zero. With the 
usual assumption of cold dark matter (w m = 0) and Q > 0, which implies that energy is 
transferred from dark energy to dark matter, some kind of exotic dark matter with a negative 
EoS is driven, assuming, of course, that we are in an expanding universe H > 0. For warm 
dark matter we would have, depending on the type of interaction considered and the strength 
of the coupling constant appearing in Q, a possible change of the sign on the effective EoS 
during the cosmic evolution. For the dark energy, the effective EoS is given by 



Q_ 

3H Pe 



WA eff = W e + nTT — ■ (2.6) 



For Q > the above equation indicates, even for w e = — 1 (cosmological constant), that 
the effective Dark Energy (DE) fluid will behave as a quintessence field. Then, an effective 
phantom behavior can only be obtained if w e < — 1. 

2.2 The behavior of the coincidence parameter parameter r 

In order to address the coincidence problem, in terms of the dynamics of the parameter 
r = Pm/Pe, w e will make in what follows a similar analysis along the line describes in ref. [12]. 
The dynamics of the parameter r is given by 

Pm_Pe\ (2?) 
Pm Pe) 

where the dot indicates derivative with respect to the cosmic time. Changing the time 
derivatives by a derivative with respect to In a 3 which will be denoted by a prime, i.e. p = 
p'3H, eqs (2.2) and (2.3) go into 

— = -(1 + ^.) + ^-, ^ = -(l + w e )--^-. (2.8) 

Pm oHp m p e SHpe 
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In terms of total density p = p e + p m , we obtain 

w m r + w e 



P 



1 + 



r + 1 



P- 



For the coincidence parameter r the evolution equation reads 

Q (1 + r) 2 " 



3H rp 



(2.9) 



(2.10) 



Before specifying any particular type of interaction, we will discuss some general of 
these equations. The critical point of Eq. (2.9) is given by 



(1 + w e ) 

W m + 1 



(2.11) 



If w m > (warm dark matter) and since r must be positive, it follows that w e < — 1, 
which corresponds to a phantom DE. Using the condition r 1 = in Eq.(2.10) and the value 
of r c one finds 



(w m - W e ) 



Q 



W m + 1 3H(1+W e ) 

Depending on the value of w m three cases can be analyze. 

2.3 Case w m = 

In this case, Eq. (2.12) simplifies to 



(2.12) 



UK 



Pc 



Q 



w m + 1 3H 



(2.13) 



This case was analyzed in [12], concluding that a positive stationary energy density 
p c requires Q > since w e < — 1. This means that independent on the interaction type, 
the existence of critical points require a positive exchange from dark energy to dark matter 
(DM). 



2.4 Case w m > (warm dark matter) 

For this case the expression for p c becomes 

W e ~ W m 



Pc 



Q 



w m + l 3H(l + w e )' 



(2.14) 



The condition w e < — 1 leads in this case to same result found in 2.3 for the sign of Q. 
This result also holds for —1 < w m < 0. 

We find here the general condition to have accelerated expansion in terms of the energy 
densities of the dark components and their EoS. Differentiating (2.1) with respect to t and 
substituting for p m and p e gives the auxiliary equation 



2H = -(1 + W m )p m - (1 + W e )p e - 



(2.15) 



The acceleration is given by the relation d = a(H + H 2 ). From (2.1) and (2.15), we 



have 



--(Pm(l + 3w m ) + (1 + 3w e )Pe)- 
D 



(2.16) 
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The condition a > leads to the condition 

l + 3w m 

P e > ~TT^ — 2 - 17 
1 + 3u; e 

Since we shall consider w e < — 1 this condition represent a right line with positive slope 
1 + 3w; m /3|w; e | — 1 in the plane (p m , p e ). 

3 Evolution equations: Fixed points and stability analysis 

Introducing the interaction term Q = Xp^pe in Eqs. (2.2) and (2.3) yields 

p m = -3(1 + w m )H Pm + Ap>f (3.1) 
Pe = -3(1 + w e )Hp e -\p^pP. (3.2) 

The time evolution of the dark matter and dark energy densities is given by the highly 
non-linear coupled differential equations (3.1) and (3.2). We rewrite them as follows: 

Pi= fi(pm,Pe), (3-3) 
where i = m,e. The functions f m and f e are defined by the following expressions 

fm(pm, Pe) = ~VS (1 + W m ) p m (p m + p e ) l/2 + Xp^Pe (3-4) 
fe(Pm, Pe) = -V3 (1 + We) Pe (Pm + Pef 1 " 2 ~ ^PmPe- ( 3 -5) 

From numerical results we expect that the above equations have some non-trivial fixed- 
points (p m , Pe), which we want to study analytically including their stability properties. In 
spite of the non-linearities and according to ref. [16], it is still possible to analyze the stability 
of fixed points by using the linearized piece of the original differential equations 

P\ = a-ijPj + Ri(Pm,Pe) (3.6) 
j 



provided the inequality 



\Ri(Pm,Pe)\<N[Y jP U (3.7) 



is fulfilled in a neighbor region around the fixed points. 

To achieve this goal we expand both the dark matter and dark energy densities around 
their fixed point values (p m , p e ) as follows: 



Pm = {Pm + PPmf = Pm + <*PPm + O (p? p m j pm) 
f4 = (p e + ppef = + P»Pe + O {p 2 pl/pe) , 

it also holds 

(Pm + Pe) 1 ' 2 = (pm + Pe^) ^ + ~ " (Pm + Pe) + O {p?) . 

2 + Pe P ) 
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Inserting the above perturbative expressions into the differential equation system one 
obtains up to first order in p: 



(3. 



or explicitly 



p e = -\/3(l +U e ) 



Pm. {Pm + P 



Pm (Pm + Pe 



a/2 



2 (Pm + Pe 



Pe {Pm + Pe) , ,_ , _ a/2 



2 {p m + Pe 



The tree level values are implicitly defined by the relations 



VS (1 + W m ) p m (p m + p e ) 1/2 = Xp^ p/ 
VS (1 + W e ) p e (p m + yO e ) 1/2 = -\p£P, 



*7i P 
e i 



or equivalently, 



l + 0J n 
l+UJe 



{-If 



A (l+a; m )' 3 - 1 (l+^ e ) 



1/2-/3 



(w e - w m ) 1/2 



(3/2-Q-/3)" 1 



/5Pe/5m + a/5m./3 e /3 (3.9) 

f3p e p^ + ap m p/] . (3.10) 



(3.11) 
(3.12) 



(3.13) 
(3.14) 



Now we are prepared to analyze the different numerical results obtained from the direct 
numerical solution of the system of eqs. (3.3). The numerical solutions were obtained by 
using a very accurate numerical method of adaptive step-size algorithm called Bulirsch-Stoer 
method, which will be explained in the next section. 

We will use the numerical values cj m = and uj e = —1-1, which are of physical interest 
as we will discuss it in the next section. For these particular values, the above equations have 
a fixed point given by p e = 10p m and 



A 10^ \ 

1 + Urn V33J 



(3/2-a-/?)- 1 



(3.15) 



These fixed points are displayed in figures 3 and 4, and their loci agree remarkably well 
with the corresponding ones of the numerical results. 

In particular, for a = 0.9, (3 = 1.0, uj m = and A = 1, eq. (3.15) leads to the relation 
p m 0.25, which is in agreement with the corresponding locus shown in figure 4. For the 
next attractor, a = 0.8, j3 = 1.0, u m = and A = 1, eq (3.15) gives p m 0.16, which again 
perfectly agreed with the numerical result shown in figure 4. 
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Now we comeback to the linearized system of eqs. (3.9) and (3.10). We rewrite it 
explicitly as the homogeneous system of differential equations 




+ Xapf p m + 



(3.16) 



WPm Pe + 



(3.17) 



The stability of the fixed points of the above equations depend on the eigenvalues given 
by characteristic equation associated to the system: 



k 2 - (an + 022) k + (ana22 - ai2«i2) = 0. 



(3.18) 



As it is well known, depending on the roots of eq. (3.18), the trajectories around the 
fixed-point (p m = 0, p e = 0) will be stable or unstable, (see for instance [16]). 

For example, if all roots of the characteristic eq. (3.18) have negative real parts, then 
the trivial solution (p m p e ) T = (0 0) T of the linearized system and of the non-linear system 
(3.3) is asymptotically stable. On the other side, if at least one of the roots of eq. (3.18) has 
a positive real part then both systems have an unstable fixed point at (0 0) T . 

For the interacting model described by eq. (3.3), there are five physical parameters, 
w m , w e , A, a and /3, which should be chosen according to both, physical stability properties 
on one side, and compatibility with observational data on the other side. The solutions of eq. 
(3.18) were numerically evaluated for different regions of the parameter space, and found the 
interesting physical region defined by the inequalities: < u) m < 1/3, — 2 — uj m < u e < — 1, 
A > 1, 0.0155 < a < 0.222 and 0.59 < a < 1.02 for uj m = and /3 > 0.8. From it, we will 
considered the subregion < uj m < 1/3, uj e = —1.1, A = 1, a = 0.9 and f3 = 1. In particular, 
in these regions the condition 3.7 holds, which guaranties that the linear analysis of stability 
also apply to the non-linear differential equation considered. 



4 Numerical analysis 

In this section we present numerical results obtained by using the Bulirsch-Stoer method to 
solve the non-linear coupled system of eqs. (3.3). The Bulirsch-Stoer method uses an adap- 
tive step-size control parameter, which ensures extremely high accuracy with comparatively 
little extra computational effort. In the past, this method has proven to be very accurate 
for solving non-linear differential equations [17]. In addition, we have computed the fixed 
point trajectories for different values of the five physical parameters w m , w e , A, a and f3. 
The corresponding trajectories within the stability region discussed in the previous section, 
will be shown in the figures bellow. 

The evolution of matter an energy densities p m and p e depends critically on the value of 
the parameter w e . In particular, for w e > —1, the system exhibits a smooth evolution of the 
densities towards the fixed point. On the other side, if w e < —1 the system shows periodic 
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orbits around the fixed point. In this case, a slight variation of the exponents a and (3 leads 
to spiral orbits as it is shown in figures 3 and 4. 

4.1 Symmetry for w e and w m 

We will now consider the particular values a = /3 = 1, which has been claimed to give the 
best fit to observations [9] for an interaction term of the form Xp^Pe- As explained above, 
we have used the Bulirsch-Stoer method to solve numerically 3.3. 

For A = 1, < w m < 1/3 and —1.0 < w e < —1.1 the trajectories of the matter and 
energy densities are displayed in the figures below. For the fixed value w m = and w e 
within the interval [—1.1, —1-0], different trajectories are shown in figure 1, which have the 
remarkable feature of having two intersection points. These points are characterized by the 
same value of p m = 0.4. For the fixed value w e = —1.1 and w m within the interval [0, 1/3], 
different trajectories are shown in figure 2. In this case the intersection points arise at the 
same value of p e = 3.5. 

In both cases the fixed points are of center point type, with the exception of w m = 
and w e = —1, which corresponds to a crossover value (see figure 1). 




Figure 1. (Color on-line) For the parameters a = j3 = A = 1, w m = and w e = [—1.0, — 1, 1] 
different evolutions of the densities are shown starting from the same initial values. Two intersection 
points can be identify for the same value of p m when the range of w e is swept. The dots represent 
the fixed points given by eq. (3.13) and (3.14). 

4.2 Spiral trajectories 

From the previous subsection, closed trajectories are obtained only for the region defined by 
w e < —1. Depending on the values of a and /3, converging or diverging spirals are obtained. 
In fact, a < 1 and /3 > 1 lead to convergent spirals and a > 1 and f3 < 1 lead to divergent 
spirals (see figures 3 and 4). 
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Figure 2. (Color on-line) For the parameters a = (3 = A = 1, w e = — 1.1 and w m = [0, — 1/3] 
different evolutions of the densities are shown starting from the same initial values. Two intersection 
points can be identify for the same value of p e when the range of w m is swept. The dots represent 
the fixed points given by eq. (3.13) and (3.14). 



5 Observational Constraints 
5.1 The Hubble parameter 

In a spatially flat FRW Universe, the Hubble constraint and the conservation equations for 
the matter and dark energy fluids are given as 

H 2 = ^f(p m + p e ), (5.1) 

p m + 3H( Pm + p m ) = Ap>f , (5.2) 

Pe + 3iT(p e +Pe) = -Ap>f, (5.3) 

where H is the Hubble parameter and A is a constant to quantify the strength of the inter- 
action between the matter with the dark energy. These equations can be written in terms of 
the scale factor a as 

r- -p m {l + w m ) = ——, (5.4a) 

da a all 

^ + - Pc (l + Wc ) = -^^. (5.4b) 
da a an 

We define the dimensionless parameter A related with A as 

_ A H 

(pi-J-Hp ^)^ 1 j 
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Q. 




Figure 3. (Color on-line) For initial values p m = 0.4 and p e = 3.5 and for the parameters a = f3 — 
A = 1, w m = and w e = —1.1, a close orbit is shown -blue line- around its fixed point -blue dot-. 
Starting for the same initial values but changing the value of the exponents beta, spiral trajectories 
are shown for the evolution of the densities -red and green lines-, for /3 = 0.9 the evolution moves 
away from the fixed points while for /3 = 1.1 the evolution is towards the fixed point . 



where p® lit = 3Hq / (8irG) is the critical density evaluated today and Hq is the Hubble con- 
stant. We define also the dimensionless parameter density fij = Pi/p% lt with i = m,e. 
Using these definitions, the Friedmann constraint equation (5.1) can be expressed as H = 
HoV&m + &c, and the conservation eqs. (5.4) become 



+ -Ml + w m ) = A — ^L=, (5.6a) 
^£ + £^(1 + Wc) = _a ^ . (5.6b) 

Using the relation between the scale factor and the redshift z given by a = 1/(1 + z) we 
express eqs. (5.6) in terms of the redshift as 



dCl m 
dz 

dCl c 
dz 



1 

l + z 
1 

l + z 



i 



A- 



3(1 + w e )tt e + A 



(5.7a) 
(5.7b) 
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Figure 4. (Color on-line) Phase diagrams for p m and p e (starting from the same initial conditions) 
are given by the color lines for different values of a, and for A = f3 = 1, w m = and w e = —1.1. They 
were obtain by solving numerically eqs. (3.4) and (3.5). For this figure and all the following ones, the 
black dashed line represent the fixed-point trajectory given in eqs. (3.13) and (3.14). The full black 
line is the zero-acceleration line given by the expression (2.15) with a = 0, the region above this line 
correspond to an accelerating universe while the region below this line correspond to a deceleration 
universe. 



We solve numerically this ordinary differential equation system (ODEs) for the functions 
f2 m (z) an d £l e (z), with the initial conditions £l m (z = 0) = O m o = 0.274, and tt c (z = 0) = 
n c0 = 0.726. 

The dimensionless Hubble parameter E = H/Hq becomes 



E(z, a, /3, A, w m ,w e ) = y Cl m (z) + (l e (z), (5.8) 
where £l m (z) and Cl c (z) are given by the solution of the ODEs (5.7). 
5.2 Cosmological probes 

We constrain the values of the free parameters (a, /3, X,w m ,w e ) using the cosmological ob- 
servations that measure the expansion history of the Universe explained in the following 
sections. We compute the best estimated value for each parameter through a minimizing 
process of a x 2 function defined below, and calculate their probability density functions (or 
likelihood function). 

5.2.1 Type la Supernovae 

We use the type la supernovae (SNe la) of the "Union2.1" data set (2012) from the Supernova 
Cosmology Project (SCP) composed of 580 SNe la [18]. The luminosity distance in a 
spatially flat FRW Universe is defined as 
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0,0 0,1 0,2 0,3 0,4 0,5 



Figure 5. (Color on-line) Phase diagrams for the parameters a = 0.8, A = /3 = 1, w m = 0.1 and 
w e = —1.1. Spiral evolution of the densities towards the fixed points is shown. 

Assumptions 



Best estimate 


Xd.o.f. 


x 2 


A 


a 


P 


w m 


W e 


A = 1.497 ±0.24 


1.30 


773.62 




0.9 


1 





-1 


A = 2.257 ±0.18 


1.29 


762.78 




0.9 


1 


0.2 


-1.1 


a = 0.379 + [j'^ 
a = 0.3Q±8:* 
f3 = 1.171 ±0.8 


1.29 
1.23 
1.22 


767.46 
732.96 
722.89 


1 
1 
1 


0.9 


1 
1 




0.1 



-1 
-1.1 
-1.1 


Wm = -0.062 ± 0.028 


1.21 


718.37 


1 


1 


1 




-1.1 


w c = -1.288 ±0.03 


1.14 


684.48 


1 


1 


1 







w e = -1.416 ±0.03 


1.10 


658.23 


1 


1 


1 


0.1 





Table 1. Best estimated values for (A, <x, fi, w m ,w e ) (left column) assuming some particular values 
indicates at the last five columns. The best estimates where computed using the joint SNe+H(z) data 
sets. See section 5.2 for details. The likelihood functions for each estimated parameter are shown in 
figures 7 to 14. as well as the stable regions of the critical points of the dynamical system analysis. 
From the value of \d f seen as measure of the goodness-of-fit to data we find that the best fit to 
the cosmological observations corresponds to the bottom row of the table where Xd o f =1-10, i.e., a 
model of a phantom dark energy (w e = —1.4) interacting with a warm dark matter (w m = 0.1), where 
the energy transfer comes from the dark energy to the dark matter (A = 1) and with an interacting 
term with the simple form "p m p c "- The errors correspond to ler (68.3%) of confidence level. 



, v c(l + z) r dz' , , 

^> = V1e(P) (5 ' 9) 

where E(z) corresponds to the expression (5.8), and "c" corresponds to the speed of light in 
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Figure 6. (Color on-line) Evolution for £l m (z) (blue) and Cl e (z) (red), with the initial conditions 
&m(z = 0) = O m o = 0.274 and Cl c (z = 0) = O e o = 0.726 (present time z = 0). The parameters are 
a = 0.9, jS = 1 and w e = —1.1. Full lines are for w m = 0.0, dashed lines arc for w m =0.1, doted lines 
are for w m = 0.2 and dash-doted lines are for w m = 1/3. The dependence of Cl (z) (red) remains 
constant with w m compare to f2 m (z) (blue). 



units of km/sec. The theoretical distance moduli //* for the k-th supernova at a distance Zk 
is given by 



//(z) = 51o£ 



d L (z) 
Mpc 



So, the x 2 function for the SNe is defined as 



XSNe(«> /5> A > '"m, W c ) = ^ 



+ 25 



(5.10) 



(5.11) 



where {if. is the observed distance moduli of the k-th supernova, with a standard deviation 
of <7fc in its measurement, and n = 580. 



5.2.2 Hubble expansion rate 

For the Hubble parameter we use 13 available data, 11 comes from the table 2 of Stern et 
al. (2010) [19] and the 2 following data from Gaztanaga et al. 2010 [20]: H(z = 0.24) = 
79.69 ±2.32 and H(z = 0.43) = 86.45 ±3.27 km/s/Mpc. For the present value of the Hubble 
parameter we take that reported by Riess et al 2011 [21] H(z = 0) = Hq = 73.8 ± 2.4 
km/s/Mpc. The x 2 function is defined as 
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Likelihood function for A 




A 



Figure 7. (Color on-line). Normalized probability density function pdf (or likelihood function) for 
the interaction coefficient A assuming the values of (a = 0.9, /? = 1, w m = 0, w c = —1) and computed 
using the SNe (blue), H(z) (red) and the joint SNe+H(z) (black) data sets. The best estimated value 
for A (the central value of the Gaussian) using the SNe+H(z) data set is A = 1.497 ± 0.24, with a 
Xd o f — 1-3- See table 1. So, it is found a non- vanishing value for the interaction coefficient, i.e., the 
interacting model is favored by the SNe+H(z) observations. This model and the positive value of A 
corresponds to an energy transfer from the cosmological constant (w e = —1) to the cold dark matter 
(w m = 0). 



Likelihood function for A 
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Figure 8. (Color on-line). Normalized likelihood function for the interaction coefficient A assuming 
the values of (a = 0.9, (3 = 1, w m = 0.2, w c = —1.1) and computed using the SNe (blue), H(z) (red) 
and the joint SNe+H(z) (black) data sets. The best estimated value for A using the SNe+TJ (z) data 
set is A = 2.25 ± 0.18, with a \\ f = 1.29. See table 1. It is found a non-vanishing value for the 
interaction coefficient. This model corresponds to a phantom dark energy (w c = —1.1) interacting 
with a warm dark matter (w m = 0.2). From the comparison of the value of \\ f = 1-29 of this 
case with that of figure 7, where Xd o f = 1-3, we find that cosmological observations slightly prefer 
the interacting 'phantom dark energy to warm dark matter model instead of a cosmological constant 
interacting with a cold dark matter. 
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H{zj) - H \ 

OH 



obs 



(5.12) 
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Likelihood function for a 




Figure 9. (Color on-line). Normalized likelihood function for the power a assuming the values 
of ((3 — 1, A = 1, w m = 0, w c = —1) and computed using the SNe (blue), H(z) (red) and the 
joint SNe+H(z) (black) data sets. The best estimated value for a using the SNe+H (z) data set 
is a = 0.37lg'ig, with a x\ o f = 1-29. Sec table 1. This case corresponds to a cold dark matter 
interacting with a cosmological constant. 

where H{zi) is the theoretical value predicted by the model and Hf hs is the observed value. 
6 Discussion and Conclusions 

In order to shed some light on the coincidence problem we have explored a cosmological model 
composed by a dark matter fluid interacting with a dark energy fluid. Motivated for very 
recent investigations we have considered a warm dark matter. Since non-linear interactions 
represent a more physical plausible scenario for interacting fluid we studied an interaction 
which is given by the term Xp^pc. We have found a general result which indicates that the 
positive critical points of the coincidence parameter r = p m / ' p e exist if w e < —1 independently 
of the interaction chosen and the particular EoS for the dark matter. We have considered 
from the beginning that the energy is transferred from dark energy to dark matter (A > 0). 

We performed an analytical analysis of the non-linear and coupled equations correspond- 
ing to the continuity equations for the dark matter and dark energy fluids. In particular we 
found the fixed points and their stability properties. 

Using a high precision numerical method we solved these equations and were able not 
only to confirm with high accuracy the analytical results but also to extend the solutions 
beyond the validity regions of the analytical analysis. 

The combined method described above allowed us to compute, in the densities space 
(Pm,Pe), the behavior of the fixed points in terms of the parameters \,a,/3,w m and w c . 
Closed orbits were found for w e < —1.1 and a = (3 = 1. If a or /3 or both of them are 
different from the unity, this closed orbits transform into spiral trajectories, evolving towards 
the origin for a < 1 and (3 = 1, and away for the fixed point a = 1 and /3 > 1. This analysis 
allowed to constrain the parameters in order to have physically reasonable scenarios, that is 
accelerated expansion in the late time phase of the cosmic evolution and far future evolution 
with finite DM and DE densities, which corresponds to spiral trajectories propagating from 
the fixed points. 

The knowledge of these stability regions was used to fit with remarkable goodness the 
observational data to the physical model proposed in this paper. 
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Likelihood function for a 
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Figure 10. (Color on-line). The upper panel shows the normalized likelihood function for the power 
a assuming the values of (/? = 1, A = 1, w m = 0.1, w c = —1-1) and computed using the SNe (blue), 
H(z) (red) and the joint SNe+H(z) (black) data sets. The best estimated value for a using the 
SNe+£f (z) data set is a = 0.30^'^, with a Xd .o .f = 1-23. See table 1. This case corresponds to a 
phantom dark energy interacting with a warm dark matter. The lower panel shows the values for a 
corresponding to stable solutions for k(a) [see eq. (3.18)], i.e., the regions where the solutions k of 
the characteristic equation (3.18) arc negatives in function of a and the given values of (/?, A, w m , w c ). 
The continuous and dashed plots correspond to the positive and negative branches of the solution of 
the quadratic equation (3.18) respectively. 



Within this stability regions, the parameters (A, a, (3, w m , w e ) were estimated through 
the cosmological observations of the Union 2.1 type la supernovae data set and the Hubble 
parameter H(z) at different redshift, that measure the late time expansion history of the 
Universe. We estimated u) m and w e to explore the nature of these fluids in the context of this 
model. To compute the estimations of each of the five parameters, we assumed particular 
values for the other four. A summary of these results and assumptions are shown on table 1 
and figures 7 to 14. 

For the interacting coefficient A, we defined a dimensionless A for convenience. Using the 
cosmological observations it is found a non-vanishing positive value for A with about 99.9% 
(4<T ) of confidence level, indicating that the data prefer the interacting model instead of the 
non-interacting one, and that the energy transfer is from the dark energy (as cosmological 
constant) to the cold dark matter, for the assumed values of (a = 0.9, j3 = 1). The same 
results are obtained when is assumed a phantom dark energy interacting with a warm dark 
matter, however, according to the values of x 2 m both case, the "phantom dark energy - 
warm dark matter" model fits better the data than the "cosmological constant - cold dark 



- 16 - 



Likelihood function for /? 



1.0 

0.8 
§ 0.6 











: SNe / 




NSNe+H(z) 


: 

Assuming: " 






\ / H(z) 


ff =0.9, u',„=0, • 








w e =-l.l, 






\y 
/\ 
/ \ 


A=l 

— 
- 



-2 2 4 6 8 10 




Figure 11. (Color on-line). The upper panel shows the normalized likelihood function for the power 
/? assuming the values of (a = 0.9, A = 1, w m = 0, w c = —1.1) and computed using the SNe (blue), 
H(z) (red) and the joint SNe+i/(z) (black) data sets. The best estimated value for j3 using the 
SNe+H(z) data set is (3 = 1.171 ± 0.8, with axj o[ = 1-22. See table 1. This case corresponds to 
a phantom dark energy interacting with a cold dark matter. The lower panel shows the values for (3 
corresponding to stable solutions for k(a) [see eq. (3.18)], where we have assumed the same values of 
(a, A, w m , w c ) as above. The continuous and dashed plots correspond to the positive and negative 
branches of the solution of the quadratic equation (3.18) respectively. The stable region for both 
branches are quite particular and narrow, however, the likelihood function for j3 using SNe+ff (z) lies 
in this stable region. 



matter" model. 

For the power a of the dark matter density, we found that a > with about 99% 
of confidence level using the SNe+-ff (z) dataset, indicating an interacting term of the form 
p e Pmi f° r the assumed values of (A = l,/3 = 1). This result is the same when we assume 
(«j m = —0, w e = —1) or (w m = —0.1, w c = —1.1), however, the best fit to data is for the case 
of a phantom dark energy interacting with a warm dark matter. 

For the power (3 of the dark energy density, we found that (3 > with about 68.3% (la) 
of confidence level using the SNe+H(z) dataset, suggesting an interacting term of the form 
A/3cPm> f° r the assumed values of (A = I, a = 0.9, w m = 0,w c = 1.1). This case corresponds 
to a phantom dark energy interacting with a cold dark matter. 

For the barotropic index w m of the dark matter EoS, we found that the best estimated 
value is negative (w m = —0.06 ± 0.03), when using the SNe+H(z) and assuming (a = 1, j3 = 
1, A = 1, w e = —1.1). When we impose the restriction w m > we simply obtain w m = 0.0. 
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Figure 12. (Color on-line). The upper left panel shows the normalized likelihood function for the 
barotropic index u? m for the dark matter, assuming the values of (a = 1, /3 = 1, A = 1, w e = —1-1) and 
the restriction w m > 0. The upper right panel shows the likelihood function for w m computed without 
assume any restriction on the value of w m . The lower panel shows the values for w m corresponding 
to stable solutions for k(w m ) [see eq. (3.18)] where we have assumed the same values of (a, j3, A, w c ) 
as above. The likelihood function of SNc +H(z) lies in the the stable region for both branches. 



And for the barotropic index w c of the dark energy EoS, we found that it is constrained 
with 68.3% (la) of confidence level to —1.32 < w c < —1.26 and —1.45 < w c < —1.39, when 
it is assumed a cold and warm dark matter respectively, using the SNe+i? (z) dataset, that 
is clearly in the phantom regime for the dark energy. 

On the other hand, from the value of Xh f seen as a measure of the goodness-of-fit to 
data, we find that from the cases studied in the present work, the best fit to the cosmological 
observations is obtained for the case with Xdof = 1-10 ( see t arj le 1) that corresponds to a 
phantom dark energy (w e = —1.4) interacting with a warm dark matter (w m = 0.1), where 
the energy transfer comes from the dark energy to the dark matter (A = 1) and with an 
interacting term with the simple form "/> m p e " • 

In summary, the non linear interaction chosen in this work leads to plausible scenarios 
that can alleviate the coincidence problem. The stable fixed points represent universes which 
end in a dark sector with non zero and finite energy densities in both fluids, despite the 
phantom behavior of the dark energy fluid. Finally, an interacting dark sector with a phantom 
like dark energy and warm dark matter fluid is well supported from cosmological observations. 
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Figure 13. (Color on-line). The upper panel shows the normalized likelihood function for the 
barotropic index w e for the dark energy, assuming the values of (a = 1, /3 = 1, A = 1, w m = 0). The 
best estimated value for w D using the SNe+if (z) data set is w c = —1.28 ± 0.03 with a x\ o f = 1-14. 
This case corresponds to a phantom dark energy fluid interacting with the cold dark matter. The lower 
panel shows the values for w e corresponding to stable solutions for k(w e ), where we have assumed the 
same values of (a, /3, A, w m ) as above. The likelihood function of SNe +H(z) lies in the the stable 
region for both branches. 
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